Valence bond entanglement entropy of frustrated spin chains 
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We extend the definition of the recently introduced valence bond entanglement entropy to arbi- 
trary SU(2) wave functions of 5* = 1/2 spin systems. Thanks to a reformulation of this entanglement 
measure in terms of a projection, we are able to compute it with various numerical techniques for 
frustrated spin models. We provide extensive numerical data for the one-dimensional Ji — J2 spin 
chain where we are able to locate the quantum phase transition by using the scaling of this en- 
tropy with the block size. We also systematically compare with the scaling of the von Neumann 
entanglement entropy. We finally underline that the valence-bond entropy definition does depend 
on the choice of bipartition so that, for frustrated models, a "good" bipartition should be chosen, 
for instance according to the Marshall sign. 
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I. INTRODUCTION 

Entanglement is a fundamental notion of quantum me- 
chanics, that has over the recent years gained popularity 
as a way to provide new insights in the quantum many- 
body problem. From the condensed matter point of view, 
one of the most interesting promises of the study of en- 
tanglement properties is the possibility to automatically 
detect the nature of quantum phases and of quantum 
phase transitions. In this approach, there is no need to 
provide a priori physical information or input, such as 
the specification of an order parameter. The detection 
can occur through the study of the scaling (with system 
size) of various entanglement estimators. For instance, 
the scaling of the von Neumann entanglement entropy 
for one-dimensional systems is different for critical and 
gapped systems - allowing their distinction. For a recent 
review of various properties of entanglement entropy in 
condenser matter, see Ref. T, 

To quantify the entanglement between two parts Q and 
Clofdi quantum system described by a wave- function |^), 
one usually invokes the von Neumann entanglement en- 
tropy (vN EE) defined as: 

S^''{n) = -Tin (pnHpn)). 

where pQ = Trj^|^)(^| is the reduced density matrix ob- 
tained by tracing out the degrees of freedom in Q. Con- 
sidering that 1] is a connected region of space (such as 
a block of sites in a lattice model), on general grounds 
one expects that S^'^ scales not as the volume of 17, but 
rather as the interface between Q and O. This "area 
law"^ is due to the fact that the value of S^'^ is indepen- 
dent of whether degrees of freedom in 17 or 17 have been 
first traced out: 5^^ (17) = S^'^iCl) and therefore the size 
dependence must come from the boundary between the 
two parts of the system. 

This general area law has been shown to be fulfilled for 
many physical wave- functions. However it is also known 
to be violated for several examples, where in most cases 



some type of long-range correlations develop between de- 
grees of freedom in Q and Cl. The most documented situ- 
ation is the case of one-dimensional (Id) quantum critical 
systems where several important results can be derived 
from Conformal Field Theory (CFT). For Id quantum 
critical wave- functions displaying conformal invariance, 
it was showiP that S^'^ = | ln(x) where x is the length 
of the block of sites ft. Here c is the central charge of the 
corresponding CFT, and periodic boundary conditions 
are assumed. The Id area law (where S^'^ saturates to 
a constant for large x) is fulfilled on the other hand as 
soon as the correlation length is finite. 

Even though one expects the area law to be valid in 
most cases, the situation of corrections to the area law is 
still unclear in higher dimensions. CFT predictions are 
valid only in a few isolated situations^ ^ . Some exact re- 
sults are available for a few specific models: for instance, 
the ground-state of free bosonic models fulfills strictly 
the area lav;^ ^ whereas free fermions can display multi- 
plicative log corrections^^. On the other hand, calcula- 
tions for interacting models become rapidly untractable. 
Numerical simulations are also not eased by the higher 
dimensionality. Exact diagonalization techniques have 
access to EE, but they are limited to very small samples 
sizes which do not allow a test of the area law, and de- 
viations thereof, in high dimension. EE comes for free 
within the DMRG methocr , but it is limited to Id and 
quasi-ld systems. Stochastic methods, such as Quantum 
Monte Carlo (QMC), have no problem with simulations 
of systems in higher dimensions, but unfortunately the 
vN EE is a quantity that is extremely complex to mea- 
sure within Monte Carlo methods (see however recent 
progresses'^^ ^^). 

Alternatively, it is possible to define for certain quan- 
tum spin systems, a different quantifier of entanglement 
through the use of the Valence Bond (VB) representation. 
The key idea is that for two quantum spins 1/2 at sites i 
and j, the singlet state (or VB) |^) = -^{\tiij) - I Utj)) 
is maximally entangled. It can therefore be used as a nat- 



ural unit of entanglement (in the quantum information 
community, the entanglement is often measured in units 
of singlets). Consider now a VB state where an even 
number N of spins 1/2 are coupled pairwise in singlets, 
and divide the spins in two arbitrary sets Q and Q. It is 
simple to see that the von Neumann entropy S^'^ is equal 
to the number of VBs shared between Q and Cl (i.e. where 
one of the two spins is located in Q and the other in Cl) , 
times the constant ln(2). This constant is just the von 
Neumann entropy of a single spin in a VB. In other words, 
every singlet that crosses the boundary between Q and Cl 
contributes ln(2) to the von Neumann entropy. The pic- 
ture is very appealing as it provides a simple geometrical 
interpretation of entanglement, a quantity which is not 
always easy to grasp intuitively. This argument of course 
is only exact for the case of VB states, which are sim- 
ple factorized states. However, it can be shown, and we 
will describe this in detail below, that the picture holds 
for all singlet states. This is an important property as 
most antiferromagnetic systems have a singlet finite size 
ground- state. 

In general, the resulting Valence Bond Entanglement 
Entropy ( VB EE)P^^ is different from the von Neumann 
entropy, except for the case of VB states where they coin- 
cide. There are however several points of interest for this 
alternative description: (i) it fulfills all desired properties 
of an entropy^, (ii) the VB EE can be easily computed 
through QMC methods in the VB basis^, offering the 
possibility to study d > 1 systems, (iii) the scaling prop- 
erties of the VB EE display several interesting features. 

Concerning the scaling (with system size) of the VB 
EE, it has been shown first numerically and then analyt- 
ically that for critical systems in Id, the VB EE scales 
logarithmically with block size S^^ ex j\n{x) (here and 
in the following the symbol oc denotes proportionality up 
to a constant). On the other hand, S^^ converges to a 
constant for gapped systems. The same scaling behav- 
ior is displayed by the vN EE. The only difference comes 
from the prefactor of the log divergence: while the numer- 
ical estimation of 7 was first reported^ to be consistent 
with c/3 (as for vN EE) , the analytical results of Ref. l15j 
indicate that the two quantities are different for Heisen- 
berg spin chains 7 = 41n(2)/7r^ 7^ 1/3, even though the 
numerical value of 7 ::^ 0.279 is very close to 1/3 (c = 1 
for the Heisenberg chain). We will comment on the nu- 
merical validation of the exact value 7 = 41n(2)/7r^ at a 
later stage of this paper. 

In two dimensions, the situation seems slightly dif- 
ferent . Fo r gapped spins systems, the VB EE was 
showrP^^ to fulfill a strict area law S^^ ex x (with x 
the linear size of the boundary between Q and Cl). The 
same scaling is expected on general grounds for S^'^ . In 
the case of the ground-state of the 2d Heisenberg model 
displaying Neel order with gapless excitations, the VB 
EE displ ays a m ultiplicative logarithmic correction to the 
area law ^^ * -^^ * -^ ^ : S^^ ex xln(x). There is no equivalent 
calculation (analytical or numerical) for the vN EE of the 
2d Heisenberg model, for the reasons described above. 



However, recent DMRG calculations^ on A^-leg ladders 
and QMC computations of the Renyi entropy of the 2d 
Heisenberg modeP^, suggest that the vN EE displays no 
such multiplicative logarithmic and that the Neel ground- 
state fulfills strictly an area law S^'^ ex x. This can be 
seen negatively as the scaling of the VB EE does not 
match the one of the vN EE, showing its limits to dis- 
cuss the adherence of the vN EE to the area law in higher 
dimensions. One should note however that the VB EE is 
able to distinguish between a gapless and a gapful state 
through its scaling -whereas the vN EE cannot-, one of 
the main original and practical motivations of studying 
entanglement in condensed-matter systems. As a side 
remark, we also note that the VB EE can be used to 
characterize shared information in the different context 
of stationary states of stochastic models^^. 

In this paper, we investigate the properties of the 
VB EE using a different approach. The original VB EE 
definition^^ is intimately related to the fact it is possible 
to consistently define^ a VB occupation number able 
to quantify the presence/absence of a SU(2) dimer on a 
given bond for any singlet state. Note that the definition 
only depends on the chosen bipartition of the lattice into 
two subsets, in spite of the VB basis overcompleteness. 
We derive in Sec. [H] an alternative but equivalent defini- 
tion of the VB occupation number which is free of any 
VB basis formulation. This allows to define the VB EE in 
the practical Sz basis, and its computation through dif- 
ferent numerical schemes (such as Exact Diagonalization 
or DMRG) than the VB QMC method used in previous 
works. Being now able to compute the VB EE for frus- 



trated systems, we study in Sec. |III| both vN and VB EE 
for the J1-J2 spin chain, using DMRG techniques. We 
discuss and compare how the scaling of both entropies 
can detect the critical phase (for small J2) and the quan- 
tum critical point that separates it from the dimerized 
gapped phase present at large J2. We also discuss the 
importance of the Marshall sign present in the ground- 
state wave function when comparing the two entropies. 
We finish with a discussion in Sec.[lVlon the usefulness of 
the approach, as well as on the further possibilities open 
by the Sz representation of VB occupation numbers. 



II. VB FREE FORMULATION OF THE VB EE 

A. VB occupation number as a projection 

Original formulation — In this paragraph we recall 
some definitions and results on VB occupation number^. 
Choosing a bipartition of the N-site lattice into two equal 
sized subsets A and S, the bipartite VB subspace is gen- 
erated by all the bipartite VB states 



\^v) 



[hj]^ 



(1) 






where [z, j] is a SU(2) dimer state and P is a dimer cov- 



ering of the system. For any bond (i, j) such as i G ^ 
and j ^ B the VB occupation number in the state |(/?x)) 
is defined as: 



nii,j){\^v)) = 



1 if {ijj) belongs to P, 
otherwise. 



(2) 



The bipartite VB manifold is overcomplete: all bipar- 
tite VB states are singlet {S = 0) states but their 
number (N/2)\ is m uch larger than the singlet sub- 
space dimensioiP^^. As a consequence, a given linear 
combination of bipartite VB states |^) = Xli^^r'l^r') 
can be rewritten in many alternative linear combina- 
tions 1^) = ^x) fJ^vl^v) with Ax) 7^ /ir). This point 
requires to reconsider the extension of Eq. Q for lin- 
ear combinations of bipartite VB states since the iden- 
tify T.v^'^^iiJ) i\^T))) = T.vl^T)n^ij) {\(pv)) is not ob- 
viously granted. It is nevertheless possible to prove that 
n(^j) is lineai^^ in \(px>) which provides an intrinsic def- 
inition of n(ij)(|^)) despite the bipartite VB manifold 
over completeness . 

Projection (VB states) — We give here an alternative 
but equivalent definition of the VB occupation number 
which (i) is explicitly independent of the way the state is 
decomposed in the overcomplete bipartite VB basis, (ii) 
is valid for any spin S and (iii) will be shown to be more 
versatile for numerical computations. The spin-S* dimer 
is defined as the two-site singlet state: 

1 +^ 

[ij]s = ^== Yl (-1)^-^1-^.,+^.) (3) 



We define the reference state 



\Rs) 



*b, -rUj . . . , O, -ta!!)/ 



(4) 



where the state is written in the ^iSz eigenstates basis 
and ordered such as A and B sites appear in alternating 
order. In particular \Ri/2) is nothing but the Neel state. 
As already noticed^^ for 6* = 1/2, the reference state has 
an equal overlap with all bipartite VB states: {Rs\^v) 
does not depend on the bipartite dimer covering V. This 
property is established by a direct evaluation from Eq. p| 
and Eq. Q showing {Rs\^v) = 1/(26' + 1)^/"^. 

For any bipartite VB state lipi)) and for any bond (i, j) 
such diS i G A and j G S we are going to show that 



Indeed, we have 



(26 + 1)^/^ 
25 



{Rs\S^S-\^v)^ (5) 



6+ 5-|i?5) = (25)|... ,-6 + 1,6-1,...) 



(6) 



If the bond (i, j) is occupied, Icpi)) = • • -^[hjjs^' • • and a 
simple inspection of Eq. isl shows that {Rs\S^S~\(pT>) = 
— (26)/(26 + 1)^/^ and hence n(^ij^{\(pT>)) as defined in 



Eq. ([5| is 1. On the other hand, \ipx>) = . . . [i, k]s 
[l^j]s... if the bond (i,jf) is unoccupied. The total 
Sz component on any occupied bond of a VB state is 
0. Thus any eigenstate of S^ + 6| (or S[ + 6|) with 
a non-zero eigenvalue is then orthogonal to |(/px)). It 

is salient from Eq. (6) that (sf^S^]s:^S-\Rs) = 



-^S^S-\Rs) and (6f 

Hence nfij\{\(pT,)) as defined in 



S^j 6/ 



S-\Rs) = -StS-\Rs). 
Eq. ([5| is in this case. 

Finally, it is easy to see that if both i and j sites are 
located in the same subset A or S, the definition Eq. ^ 
also ensures that ^(i,j) = which is always true (in- 
dependently of \(pv)) as no dimer is allowed on such a 
non-bipartite bond. 

Let us mention some of the advantages of definition 
Eq. ([5| as an alternative to Eq. ([2|. First of all, it is 
explicitly linear in \(fT>) which ensures that its extension 
to arbitrary linear combination of bipartite VB states can 
be consistently defined. As stated before, in the case of 
6 = 1/2, the subspace of bipartite VB states is a basis of 
the total singlet sector, ensuring that a dimer occupation 
number can be defined for any 6 = 1/2 singlet. This is 
not true anymore for general spin 6: bipartite spin-S VB 
states do not form a basis of spin-6 singlets. However 
it can be shown that they form a basis of the subset 
of spin-S singlets that are also SU(7V) singlet^^ (with 
A^ = 26 + 1) so that Eq. ([5| can be useful in that context. 

Projection (VB states superpositions) — Using Eq. ([5| , 
the occupation number for an arbitrary linear combina- 
tion of bipartite VB states 
as. 



1^) = ZId Ax)|(/9d) is defined 



n^^,J)m) = -^ 



(26 + 1)^/^ 
26 



{Rs\S^S-\^), (7) 



where A/" is a normalization constant. It would be tempt- 
ing to take JV = l/^/{^[^. However n(^ij^ is designed 
to measure the number of dimers (0 or 1) on bond {i,j) 
and for any VB state a given site i G Ais dimerized with 
another site on sublattice B. Hence the normalization 
condition writes, 



^n(,,,)(|*)) = l. 
iee 



This condition enforces, 
1 



Af = 



EpAx, {2S + ir/^Rs\^)' 



and 



n(i,,)m) = 



1 {Rs\S^S-\^) 
"25 (i?s|*) 



(8) 



(9) 



(10) 



Contrary to Eq. Q, this last expression is explicitly 
independent of the linear combination chosen to expand 
1^) on the overcomplete bipartite VB manifold as it only 
involves projections of |^). Since Eq. (ItI does not give 
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FIG. 1. (Color online) Three possible cases for valence bond 
configurations in the non-crossing basis for a given block B 
with an odd number of sites and neighboring sites L and R 
(see text for details). 



any prominent role to the bipartite VB basis to express 
1^), it will allow numerical computations outside the VB 
QMC scheme su ch a s with Exact Diagonalization and 
DMRG (see Sec.[lIB|. 

Note that this expression is potentially singular if \Rs) 
is orthogonal to |^). As an example, let us consider a 
one-dimensional spin- 1/2 chain with N = Ap sites (where 
p is an integer) . If we denote S the spin inversion sym- 
metry S^ — ?► —S^ and T the translation symmetry, any 
q = t: singlet state |^) will transform as 5|^) = |^) and 
T\^j) = -|^). Consequently, 5T|^) = -|^). On the 
other hand, if the bipartition ABAB ... is chosen, the 
reference state given in Eq. Q is obviously invariant un- 
der ST- As a consequence, |^) and \Rs) are orthogonal 
and Eq. (10) can not be used. 



This issue, which is a direct consequence of the normal- 
ization defined by Eq. |9|, suggests that the bipartition 
and hence \Rs) (see Eq.^ may not be chosen regardless 
of 1^). More generally, normalizing a state by the sum 
of its coefficients in an expansion like in Eq. ^ , requires 
a careful inspection of its nodal structure or equivalently 
of its Marshall sign which in turn dictates an appropri- 
ate bipartition for the reference state. We will further 
discuss this issue in Sec. IIII CI 

VB EE — Using Eq. (IlO]), the VB EE measuring the 



entanglement between 1] or 1] in a state |^) can be ex- 
pressed as. 



si'^m) 



ln2 X 



E 

(i,j) such as 



Hi,i)m). (11) 



where the spatial sums run over all possible locations of 
VBs, that is over all sites imO. and all sites j inCl. Since 
we know that n(i^j) = whenever i and j belong to the 
same subset A or S, we can restrict the summation only 
to the non-vanishing cases. 



VB EE for one- dimensional periodic systems — We 
finally make a remark on the behavior of S^^ for 
translation-invariant one-dimensional systems. When Vt 
is a linear segment of size 2n (with n integer) we find 
that 

5^^(2n) = i [5^^(2n - 1) + 5^^(2n + 1)] , (12) 

which means graphically that S^^ is made of linear seg- 
ments. Let us propose an easy graphical proof of this 
statement. Consider a VB configuration and let us com- 
pare VB EE for different blocks obtained by adding two 
extra sites R and L at each end of a (2n — l)-sites block 
B (see Fig. fl]). Since the VB EE is a well defined quan- 
tity, we can choose to work in the complete non-crossing 
basis^^ where VB do not cross according to some one- 
dimensional ordering of the sites. Since the 2n — 1 block 
has an odd number of sites, it is clear that R and L can- 
not be connected by a singlet. Then, we can consider 
all possible cases: (i) R is connected to B but not L, 
or vice- versa; (ii) neither R nor L are connected to B] 
(iii) both R and L are connected with B. These cases 
are shown in Fig. [l] and from the figure, it is straight- 
forward to check that S"^^{B) -^ S^^{B ^ R ^ L) = 
S'^^iB + i?) + S^^{B + L). As a conclusion, if periodic 
boundary conditions are used so that the entropy only 
depends on the number of sites of the block, we deduce 
Eq. ([12). 



B. Numerical computations 

From now on, we focus on the case of 5 = 1/2 systems. 

With Exact Diagonalization — We use the Lanczos 
algorithm in order to compute the ground-state of large 
Id chains^^ . We also implement lattice translations as 
well as fixing the total S^ quantum number in order to 
reduce the Hilbert space size so that we can solve systems 
up to N = 32 sites. Once the wavefunction is obtained in 
the symmetrized basis, one can easily compute its overlap 
with the reference state, which gives the denominator of 
Eq. (10). In order to compute its numerator, we need to 



apply the operator S^S, for all pairs of sites (i^j) with 
i in the selected block and j outside (let us remind that 
we can restrict ourselves to the case where i and j belong 
to different subsets). In this case, it turns out that it is 
simpler to apply this spin operator on the reference state, 
since it reduces to a swap operator for spins 1/2. Finally, 
in order to compute the VB EE, and since we are using 
translation symmetry, we have to make an average over 
all the positions of the block on the chain. 

With DMRG — The calculation of the VB entangle- 
ment entropy is straightforward, utilizing matrix product 
techniques to calculate the overlap between the ground- 
state and the reference state N = (m ... i|^), and 
the expectation value P = (t|t ... | |7^|^) where 
V = ZliGQjG^ '^t'^i + ^/^z~ h^s ^ simple representa- 
tion as a Matrix Product Operator, using the techniques 



described in Ref. ^23^, The VB entanglement entropy is 
then S^^ = - \n{2)P/N. 



III. RESULTS FOR THE J1-J2 SPIN CHAIN 
A. Model and simulation details 

We now present numerical results for the frustrated Ji- 
J2 spin chain. |^) in Eq. ( 11 ) is taken as the ground-state 
of the S = 1/2 Hamiltonian : 



H 



L 



Jl^i'^i-\-l + J2^i'^i-\-2 



(13) 



where we set Ji = 1 and will vary J2. The physics of 
this spin chain is well understood: for J2 smaller than 
the critical value J| c^ 0.241167 (Ref. ^, the system 
displays antiferromagnetic quasi-long range order, with 
algebraically decaying spin correlations. For J2 > Jf, 
the system is located in a gapped dimerized phase which 
spontaneously breaks translation symmetry. We will 
study both VB and vN EE entropies in this system in 
both phases. 

Results for J2 < 0.5 were obtained with DMRG. We 
used samples with L = 64, 128 and 192 and periodic 
boundary conditions in order to avoid dimerization ef- 
fects in the entanglement entrop}'^^ ^^, which complicate 
the finite-size analysis. Up to m = 1092 SU(2) states 
(roughly corresponding to 4000 usual U(l) states) have 
been kept for the largest samples. A long warm-up pro- 
cedure has been used, by performing between 10 and 50 
sweeps each time m was increased by 50. Convergence 
has been checked by ensuring that the energy does not 
change significantly on more than 20 sweeps for the last 
value of m. Truncation error per site and variance per site 
{H — E)'^/L were always at most 10~^^ for the largest sys- 
tems. For these periodic boundary conditions, a two-sites 
version of the DMRG algorithm has been used. We will 
essentially present results for the largest L = 192 chains, 
but will occasionally show data for smaller L when a dis- 
cussion of finite-size effects is necessary. 

Prior to calculating the scalar products with the refer- 
ence state, we use the Wigner-Eckart theorem to project 
the SU{2) groundstate to /7(1), thus giving direct access 
to the axis-dependent spin vector operators. 

Later in the paper, we will present results for J2 > 0.5 
which were obtained with ED for chains of length up to 
1/ = 32. For large values of J2, the DMRG algorithm 
has more difficulties to converge, even for small samples 
- a fact which has already been reported^. Also, some 
intrinsic difference shows up in the definition of the VB 
EE in this case due to the rapid vanishing of the Marshall 
sign in the ground-state wave function. In that situation, 
the analysis of the definition as well as meaning of the 
VB EE is different and will be discussed in Sec. IIIICI 

We finally note that the vN EE of the ground-state of 
this spin chain was studied previousl}/^^, albeit on smaller 



systems, with ED techniques. As we will see later, the 
use of large systems is necessary to locate precisely the 
quantum phase transition at J^. 



B. Results for J2 < 0.5 

We will present in this section results obtained for 
vN and VB entanglement entropies in parallel. We first 
present raw data for both entropies as a function of the 
block size x for different values of J 2 in Fig. [2j Data are 
shown only for x G [0, L/2] (we have checked that curves 
are symmetric around L/2). Both sets of curves show 
a similar behavior: on the scale of the figure, one can 
distinguish between curves which converge to a constant 
(for J2 > 0.4) and those which grow slowly but steadily 
with X. The difference between the two entropies appears 
on the former cases, where curves for different J2 appear 
more shifted for S^^ . The shift also exists for S^^ but 
is smaller (see zoom around x ~ L/2 in the inset of the 
figure) . 

One should also note the clear dimerization of both 
entropies for large J2: this is naturally expected at the 
Majumdar-Ghosh29 point J2 = 1/2 where S^^ and S'^^ 
are strictly equal to for even x and ln(2) for odd x. 
Note that this dimerization effect comes from the intrin- 
sic dimerized nature of the ground state in this region, 
and not from the boundary conditions as in Ref. 27, 

Let us concentrate now on the upper beam of curves, 
for J2 < 0.35. From conformal invariance of the ground- 
state in the critical phase, the use of the conformal block 
length x' = L/7Tsm{7Tx/L) should be useful for systems 
with periodic boundary conditions: in the critical phase, 
vN EE should scale as S^^ = c/3\n{x') -\- Ki whereas 
S^^ =-f\n{x')^K2. 

Fig. [3] displays both entropies versus x^ in a log-linear 
scale. All curves seem at first glance linear with approx- 
imately the same slope, except for J2 = 0.35 where a 
crossover to a constant regime can be identified: this is 
well visible for S^^ in the figure, but is also the case for 
S^^ when zoomed in. 

We fit the curves to a form 6'^^ = Ceff/31n(xO + Ki 
and S^^ = 7eff In(x^) -h K2 within the window x^ > 10. 
Fits are excellent and lead to Ceff (respectively 7eff) very 
close to the CFT prediction 1 (resp. 41n(2)/7r^) in the 
critical phase. The values of the fitted Ceff and 7eff are 
displayed in the left insets of Fig. [3J as obtained for the 
three different samples sizes L used in this study. The 
finite-size effects are found to be small on both quantities. 

Several remarks are in order at this stage : 

(i) The finite size dependence of the fitted values indi- 
cate that J2 = 0.35 is clearly not in the critical regime 
as we already guessed from a visual inspection of curves. 

(ii) It is quite interesting to note that both 7eff and Cgfr 
are not strictly equal to the predicted values for low J2 
(including J2 = 0) but are getting closer monotonously 
to the theoretical predictions when increasing J2. This 
effect can clearly be seen for 7eff, but also exist (even if 
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FIG. 2. (Color online) Scaling of von Neumann S^^ (top panel) and Valence Bond 5*^^ (bottom panel) entanglement entropies 
as a function of block size x, for different values of the frustrating coupling J2. Chain length is L = 192. Top inset is a zoom 
of 5'^^ for X close to L/2. 
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FIG. 3. (Color online) Scaling of von Neumann S^^ (top panel) and Valence Bond 5*^^ (bottom panel) entanglement entropies 
as a function of conformal block size x' — L/7Tsm{7Tx/L) in a log-linear scale, for different values of the frustrating coupling J2. 
Chain length is here L — 192. Top left (bottom right) inset shows the value of the coefficient of a log fit for both von Neumann 
(Valence Bond) entanglement entropies as a function of J2 , for different system sizes (see text for details) . Top right inset gives 
the value of the non-universal additive constant Ki obtained from the same fit of S^^ , as a function of J2 for L = 192. 



small) for Ceff. 

(in) The values closest to 1 and 41n(2)/7r^ are found to 
be precisely at J^ = 0.241167. The fitted Ceff is smaller 
than 1 for J2 = 0.30, and 7eff is smaller than 41n(2)/7r^ 
for J2 > 0.241167. 

The two former points lead to the following interpre- 
tation : given the existence of dangerously irrelevant op- 
erators in the critical phase (but not at J^), we expect 



that they could influence the effective value of the central 
charge and 7 as measured from a fit of EE on finite sys- 
tems (they should not in the thermodynamic limit). The 
strength of their influence decreases as one approaches 
the critical point where they vanish. This scenario sounds 
plausible for 7: indeed in the field-theoretical descrip- 
tion of the Heisenberg chain, 7 is related to the cou- 
pling constant of the free boson field whose numerical 
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FIG. 4. (Color online) Difference of the entanglement entropies to the predicted analytical values S^^ — l/31n(x^) (top panel) 
and S^^ — 7ln(x^), as a function of conformal block size x\ for different values of the frustrating coupling J2. Chain length is 
here L = 192. 



determination is known to suffer from log corrections 
due to dangerously irrelevant operators. A similar ef- 
fect has been recently predicted for the effective central 
charge in presence of marginally irrelevant operators^, 
with the prediction that Ceff < c. Raw fits of the form 
5*^^ = Ceff/31n(x') + Ki appear to give a value of Ceff 
slightly larger than 1 (in all cases less than 1%). We 
find however that the simultaneous fits of Ceff and Ki 
actually affect the determination of Ceff . A more precise 
fitting procedure^ along the lines of Ref. [31] produce val- 
ues of Ceff < 1 in agreement with Ref. 30| Details of such 
a precise estimation of Ceff are left for a future study (we 
checked that a similar analysis for 7eff does not affect the 
results displayed in Fig. |3|. 

The analysis above explains why first simulations of 
the unfrustrated Heisenberg chain (at J2 = 0) indicate 
that the scaling of the VB EE was identical to the one of 
vN EE: indeed the fitted value of 7eff ~ 0.310 is closer to 
1/3 than to 4 ln(2)/7r2 c± 0.281. We note that the transfer 
matrix estimates of 7eff also display such a small discrep- 
ancy for the Heisenberg chain^^. The numerical results 
of Ref. [15] for other spin chains not suffering from these 
log corrections appear to be in much better agreement 
with the analytical predictions, confirming this scenario. 

Actually, the vanishing of these log corrections appear 
as a way to detect on finite systems the quantum critical 
point J2 through the log scaling of 5*^^ (and possibly 
5*^^), as long as the exact values 7 = 4 ln(2)/7r^ and c = 1 
are known. If these values were not available, it would be 
more difficult to judge on the extent of the critical phase. 
Indeed from the sole quality of the fits, data at J2 = 
0.25 and J2 = 0.30 (which are theoretically located in 
the gapped phase) are compatible with a critical scaling. 
This is certainly due to the small simulation length used 
L with respect to the large correlation length close to Jf . 

Finally, we discuss the behavior of the constants Ki 



and K2. Both constants are non-universal and are a pri- 
ori not related. The fitted value of Ki is shown in the 
right inset of Fig. |3J and displays a non-monotonous be- 
havior (especially at small J2). This non- monotonous 
behavior can also directly be seen on the raw 5*^^ data 
at L/2 (inset of Fig. [2]) . On the other hand, and this 
can be noticed without a fit, the constant K2 for S^^ 
decreases monotonously with J2. 

The final transformation which summarizes these re- 
sults consists in directly subtracting the expected exact 
value from both entanglement entropies. Fig. |4] displays 
gvN _ i/3in(x') and 5'^^ - 7ln(xO versus x^ Curves 
should saturate to the constants Ki and K2 respec- 
tively, which they do (except obviously for J2 = 0.35) 
on this scale. Zooming in, one observes that all curves 
for S^'^ grow in a very smooth way, except for J2 = 0.30 
which actually decreases with x^ (for 5^^, all curves for 
J2 > 0.241167 tend to decrease when increasing block 
size x^). This is in correspondence with the fitted Cgff < 1 
for this value of J2. The flattest curves are observed for 
J2 = 0.241167 ~ J2 and J2 = 0.25 for both entropies, in 
agreement with our previous observations. 



C. Results for J2 > 0.5 

Knowing the nodes and the signs of all coefficients of a 
correlated wavefunction is a difficult task. Indeed, such 
an information could allow to design a sign-free QMC 
algorithm as well as help on building variational wave- 
functions. In this context, one can define the so-called 
Marshall-Peierls sign^ 



Sim 



^{-l)Ka,\a, 



(14) 



AB bipartition 
i f t M i ^ 1 ■ r 




1 


AABB bipartition 


0.8 


jj 1 ' 1 -^^ 


0.6 


j\l^^ 


0.4 


it 


0.2' 


I 


n 


I.I. 



1 2 

J2 



FIG. 5. (Color online) Marshall sign vs J2 for a L = 32 
chain for both the lowest ^ = and q — n eigenstates. Left : 
AB partition for the Marshall sign. Right : AABB partition. 
Ground- state results are denoted with filled symbols. 



where the sum runs over the IV^i) Sz basis states and 
with the wavefunction given by |^) = X^^ttilV^i). N^ 
counts the number of up spins on the B sublattice so 
that obviously, s depends on the choice of bipartition. 

For non-frustrated Heisenberg model, it can be shown 
that the ground-state has 5 = 1 with the natural choice of 
bipartitioiP^. Frustration will spoil this result, although 
the sign may not drop suddenly (see for instance Id or 
2d Ji — J2 model^^). On left panel of Fig. pj we present 
ED data for the frustrated Heisenberg chain. With the 
natural AB bipartition, the Marshall sign stays extremely 
close to 1 for < J2 < 0.5 but starts to deviate substan- 
tially beyond. Since the ground-state oscillates between 
having momentum and tt, we plot both values of s. 

Fig. [6] shows ED data for various entropies. In par- 
ticular, since S^^ depends on the choice of bipartition, 
one may wonder what to choose. Usually, one is guided 
by the Ising solution: for small J2/J1, it is natural to 
choose a ABAB . . . bipartition, while for large J2/J1, the 
system will behave as two decoupled Heisenberg chains 
with twice as large lattice spacing, meaning that biparti- 
tion should be of the form AABB . . . The Marshall sign 
for this AABB partition is presented in the right panel 
of Fig. [S] 

Of course, the intermediate region with maximal frus- 
tration has no preferred bipartition. Moreover, since the 
ground-state oscillates between q = {) and q = n^ we plot 
both entropies. Still, as can be seen from its definition 
Eq. (11), S^^ is only defined when the overlap between 



the ground-state and the classical Neel state is finite. 
Unfortunately, the ABAB Neel state has no projection 
in the singlet q = it sector as it is invariant by a combi- 
nation of lattice translation and spin reversal. However, 
for AABB bipartition, we can compute S^^ for both low- 
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the ground- state. 



est q = and q = tt states and it turns out that data 
are very similar (sometimes they cannot be distinguished 
on the scale of Fig.[6|, although states are quite different 
(see their Marshall sign in Fig. Isl. 

By comparing Fig. |6] and Fig. |5J we observe that when 
the Marshall sign is too small, S^^ has absolutely no 
meaning (and can even be negative). On the other hand, 
is s is large enough, or said differently, if we choose the 
bipartition that maximizes 5, then S^^ behaves much 
better and follows the same trend as 6*^^, that is both 
converge to a constant for large enough block size. In 
fact, for large J2/J1 where the Marshall sign becomes 
again close to 1 for the AABB choice, we observe that 
both entropies become more and more similar. 



IV. DISCUSSION 

In this paper, we studied the behavior of the Valence 
Bond Entanglement Entropy of the frustrated J1-J2 spin 
chain, and offered a direct comparison with the von Neu- 
mann entropy. Numerical DMRG calculations indicate 
that both entropies scale logarithmically with block size 
in the critical phase, and converge in the gapped dimer- 
ized phase of this model. The study of the VB EE has 
been made possible in this frustrated model through a 
formulation of valence bond occupation number, which 
extends this notion out of the Valence Bond basis. 



We now discuss several interests of studying entangle- 
ment in quantum spin systems through a valence bond 
measure, and point out some open issues. 

First, as based on the example of the J1-J2 model, the 
scaling of S^^ with the block size allows to differentiate 
critical from gapped phases in one dimension, similarly 
to S^^ . Moreover, the knowledge of the exact prefactoi^^ 
of the log scaling permits a relatively precise determina- 
tion of the quantum critical point at J^ with finite-size 
data (this is however due to the vanishing of log cor- 
rections at this particular point, a non- generic feature). 
Note that this knowledge can be useful as the scaling 
of S^^ might now be used to characterize uniquely the 
unknown phase of a new model. S^^ can be computed 
for the Q-states Potts model using the loop language of 
Ref. [15] and there, the prefactor of the logarithmic scal- 
ing depends on Q and is therefore indeed different for the 
different critical points encountered in the Potts model. 
This is similar to the scaling of S^^ which allows the 
determination of the central charge - the knowledge of 
which entirely determines the CFT at play for minimal 
models. 

It is possible to compute S^^ directly in the thermody- 
namic limit using the infinite-size algorithm iDMRG^^ or 
iTEBD^^. In this formulation, the number of basis states 
in the calculation controls the spectrum transfer matrix 
of the system, which gives a scaling of the correlation 
length ^ (X m'^ at criticality, where a>: is a function of the 
central charge^^. However, the effective boundary condi- 
tion for the transfer matrix, and hence the form of the 
corrections to scaling of the entropies S^^ and S^^ ^ will 
be different to the case of periodic boundary conditions, 
and we leave this analysis for a future study. Besides, 
knowing the corrections to scaling induced on a finite 
system by marginally irrelevant operators is interesting 
by itself (see for instance Ref. 30). 

In dimension higher than 1, it was already demon- 
stratecP^n^nSl that the scaling of S^^ discriminates be- 
tween gapped and gapless phases. Since numerical calcu- 
lations are possible in d > 1 with QMC VB methods^^, it 
would be of high interest to perform a systematic study 
of the scaling of S^^ in different phases of quantum spin 
models. Several questions are in order: for instance, 
do the multiplicative log corrections observed for the 2d 
Heisenberg model have a physical interpretation? Are 
prefactors of the scaling of S^^ universal within a phase 
or at a quantum critical point^^, as observed in Id? 

We note that the techniques described here for calcu- 
lating the VB entanglement entropy can easily be applied 



to higher dimensional tensor network algorithms such as 
PEPS^^ where the necessary scalar products are similarly 
easily computed. This opens the door to studies of the 
VB entanglement entropy for frustrated 2D systems^. 

Our study on frustrated systems also sheds lights on 
the importance of a good physical choice for the bipar- 
tition used in the definition of the VB EE, and its rela- 
tion to the existence of a Marshall sign rule (or a large 
Marshall sign) in the wave-function under study. When 
the Marshall sign is exactly or close to 1, the resulting 
choice of bipartition (or equivalently reference state \Rs) 
in Eq. (10)) produces a Valence Bond entanglement en- 
tropy that closely follows the scaling of the von Neumann 
entropy. This suggests another route to quantify entan- 
glement in a wave- function through its projection over a 
well-chosen (physical) reference state. 

Another interesting situation which we have not dis- 
cussed in this work is the one of strongly disordered spin 
systems, where S^^ and S^^ coincide^^ (after averaging 
over disorder). This is the case in the random singlet 
phase, where the low-energy physics is dominated by a 
single valence-bond state, as remarked in Ref. 12 and 
more recently by Tran and Bonesteel^^. As pointed out 
by these authors, the study of the fluctuations of the 
number of VBs crossing the boundary provide additional 
physical insights in this situation. 

Finally, we comment on the usefulness of the formu- 
lation Eq. (10) of the valence bond occupation number. 
It clearly points towards generalizations of S^^ for situ- 
ations not explored before, for instance for spins higher 
than 1/2, as well as for systems which lack SU(2) symme- 
try. In the latter case, the direct interpretation in terms 
of SU{2) VBs is not possible anymore and the physi- 
cal meaning of 7i(^ j) has first to be clarified. It would 
be interesting for instance to look for the relation to q- 
deformed singlets^, which are used in Ref. [15] to extend 
the VBEE to Potts models. Another high interest of 



Eq. ( 10 ) is that it allows analytical insights on the distri- 



bution of valence bonds and their correlations in a singlet 
ground-state^. 



ACKNOWLEDGMENTS 

We warmly thank O. Giraud for a remark that led to 
this work, M. Oshikawa for collaboration on related work 
and P. Calabrese and F. Becca for very useful discussions. 
We thank CALMIP for allocation of CPU time. This 
work is supported by the French ANR program ANR-08- 
JCJC-0056-01. 



^ P. Calabrese, J. Cardy, and B. Doyon, J. Phys. A 42, 

500301 (2009). 
^ M. Srednicki, Phys. Rev. Lett. 71, 666 (1993) ; J. Eisert, 

M. Cramer and M.B. Plenio, Rev. Mod. Phys. 82, 277 



(2010). 
^ C. Holzhey, F. Larsen, and F. Wilczek, Nucl.Phys.B 424, 
443 (1994); C. Vidal, J.I. Latorre, E. Rico and A. Kitaev, 
Phys. Rev. Lett. 90, 227902 (2003); J.I. Latorre, E. Rico 



10 



and G. Vidal, Quant. Inf. Comput. 4, 48 (2004); P. Cal- 

abrese and J. Cardy, J. Stat. Mech. P06002 (2004). 
^ E. Fradkin and J.E. Moore, Phys. Rev. Lett. 97, 050404 

(2006). 
^ J.-M. Stephan, S. Furukawa, G. Misguich and V. Pasquier, 

Phys. Rev. B 80, 184421 (2009). 
^ M. Cramer, J. Eisert, M. B. Plenio and J. Dreissig, Phys. 

Rev. A 73, 012309 (2006). 
^ T. Barthel, M.-C. Chung and U. Schollwock, Phys. Rev. A 

74, 022329 (2006). 
^ M.M. Wolf, Phys. Rev. Lett. 96, 010404 (2006); D. Gioev 

and L KUch, ihid 96, 100503 (2006). 
^ S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 
^° P.V. Buividovich and M.L Pohkarpov, Nucl. Phys. B 802, 

458 (2008). 
^^ M. Caragho and F. Ghozzi, JHEP 11 (2008) 076. 
^^ F. Alet, S. Capponi, N. Laflorencie and M. Mambrini, 

Phys. Rev. Lett. 99, 117204 (2007). 
^^ R.W. Chhajlany, P. Tomczak and A. Wojcik, Phys. Rev. 

Lett. 99, 167204 (2007). 
^^ A.W. Sandvik, Phys. Rev. Lett. 95, 207203 (2005); A.W. 

Sandvik and H.-G. Evertz, Phys. Rev. B 82, 024407 (2010). 
^^ J.L. Jacobsen and H. Saleur, Phys. Rev. Lett. 100, 087205 

(2008). 
^^ A.B. Kahin, L Gonzalez, M.B. Hastings and R.G. Melko, 

Phys. Rev. Lett. 103, 117203 (2009). 
^^ M.B. Hastings, L Gonzalez, A.B. Kalhn and R.G. Melko, 

Phys. Rev. Lett. 104, 157201 (2010). 
^^ F.C. Alcaraz, V. Rittenberg and G. Sierra, Phys. Rev, E 

80, 030102(R) (2009); F.C. Alcaraz and V. Rittenberg, J. 

Stat. Mech. P03024 (2010); P. Calabrese J. Stat. Mech. 

N05001 (2010) ; A.L. Owczarek, J. Stat. Mech. P12004 

(2009). 
^^ M. Mambrini, Phys. Rev. B 77, 134430 (2008). 
^° G. Rumer, E. Teller and H. Weyl : Nachr. Gott., Math- 

physik. Klasse, 499 (1932). See also H. N. Temperley and 

E. H. Lieb, Proc. Roy. Soc. Lond. A. 322, 251-280 (1971) 



and R. Saito J. Phys. Soc. Japan 59, 482-491 (1990). 
^^ K. S. D. Beach, F. Alet, M. Mambrini and S. Capponi, 

Phys. Rev. B 80, 184401 (2009). 
^^ See for instance N. Laflorencie and D. Poilblanc, Lect. 

Notes Phys. 645, 227 (2004) and references therein. 
22 LP. McCulloch, J. Stat. Mech. P10014 (2007). 
2^ R.W. Chhajlany, P. Tomczak, A. Wojcik, and J. Richter, 

Phys. Rev. A 75, 032340 (2007). 
2^ W. Marshah, Proc. R. Soc. London Ser. A 232, 48 (1955); 

E.H. Lieb and D.C. Mattis, J. Math. Phys. 3, 749 (1962). 
2^ S. Eggert, Phys. Rev. B. 54, R9612 (1996). 
2'^ N. Laflorencie, E.S. S0rensen, M.-S. Chang and L Affleck, 

Phys. Rev. Lett. 96, 100603 (2006). 
2^ S.R. White and L Affleck, Phys. Rev. B 54, 9862 (1996). 
2^ C.K. Majumdar and D.K. Ghosh, J. Math. Phys. 10, 1388 

(1969); ihidnm (1969). 
2° J. Cardy and P. Calabrese, J. Stat. Mech. P04023 (2010) . 
2^ J.C. Xavier, Phys. Rev. B 81, 224404 (2010). 
^^ P. Calabrese, private communication. 
22 E.H. Lieb, T.D. Schultz, and D.C. Mattis, Ann. Phys. 

(NY) 16, 407 (1961). 
2^ J. Richter, N.B. Ivanov, and K. Retzlaff, Europhys. Lett. 

25, 545 (1994). 
25 L P. McCulloch, l arXiv:0804.2509 l (unpublished). 
2^ G. Vidal, Phys. Rev. Lett. 98, 070201 (2007). 
2'^ F. Pollmann, S. Mukerjee, A. M. Turner and J. Moore, 

Phys. Rev. Lett. 102, 255701 (2009). 
2^ F. Verstraete and J.L Cirac, cond-mat /0407066| F. Ver- 

straete, M. M. Wolf, D. Perez-Garcia, and J. L Cirac, Phys. 

Rev. Lett. 96, 220601 (2006). 
29 V. Murg, F. Verstraete and J. L Cirac, Phys. Rev. B 79, 

195119 (2009) 
^° G. Refael and J.E. Moore, Phys. Rev. Lett. 93, 260602 

(2004). ^ 



H. Tran and N.E. Bonesteel, arXiv:0909.0038| (unpub- 
lished). 

V. Pasquier and H. Saleur, Nucl. Phys. B 330, 523 (1990). 
M. Oshikawa, D. Schwandt, and F. Alet, unpublished. 



